arXiv:1504.01082v2 [cond-mat.mtrl-sci] llJul2015 


Enhancing thermoelectric figure-of-merit by 
low-dimensional electrical transport in phonon-glass 

crystals 


Xue-Ya Xiaoxiang Kai-Lun Yao,^ Xiaoming Huang,Nuo Yang,*’^’§ 

and Jing-Tao 

School of Physics and Wuhan National High Magnetic Field Center, Huazhong University of 
Science and Technology, Wuhan 430074, P. R. China, Nano Interface Center for Energy (NICE), 
School of Energy and Power Engineering, Huazhong University of Science and Technology, 
Wuhan 430074, P R. China, School of Energy and Power Engineering, Huazhong University of 
Science and Technology, Wuhan 430074, P R. China, State Key Laboratory of Coal Combustion, 
Huazhong University of Science and Technology, Wuhan 430074, China, and Beida Information 

Research, Tianjin 300457, China 

E-mail: nuo@hust.edu.cn; jtlu(S)hust.edu.cn 

Abstract 

Low-dimensional electronic and glassy phononic transport are two important ingredients 
of highly-efficient thermoelectric material, from which two branches of the thermoelectric re¬ 
search emerge. One focuses on controlling electronic transport in the low dimension, while 
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the other on multiscale phonon engineering in the bulk. Recent work has benefited much 
from combining these two approaches, e.g., phonon engineering in low-dimensional materials. 
Here, we propose to employ the low-dimensional electronic structure in bulk phonon-glass 
crystal as an alternative way to increase the thermoelectric efficiency. Through firsf-principles 
elecfronic sfrucfure calculation and classical molecular dynamics simulafion, we show fhaf fhe 
71-71 Slacking Bis-Dilhienolhiophene molecular cryslal is a nalural candidafe for such an ap¬ 
proach. This is determined by fhe nalure of ils chemical bonding. Wifhoul any opfimizalion of 
fhe material parameter, we obfain a maximum room-lemperafure figure of meril, ZT, of 1.48 
af opfimal doping, fhus validafing our idea. 

Keywords: Thermoeleclric Effecl, Molecular Crysfals, Elecfronic Sfrucfure Calculation, Molec¬ 
ular Dynamics Simulation 


Introduction 

The thermoelectric efficiency of a material is characterized by the dimensionless figure of merit 
ZT = oS^T /(tQ + Kph), with S the Seebeck coefficient, o the electrical conductivity, the elec¬ 
tron thermal conductivity, Kph the phonon thermal conductivity, and T the temperature. The opti¬ 
mization of ZT is a highly non-trivial task, e.g., the electrical conductivity and the electron thermal 
conductivity are, in many cases, related by the Wiedemann-Franz law. 

Different strategies have been proposed to increase the ZT value.^^ In their seminal work, 
Hicks and Dresselhaus proposed to increase S utilizing the sharp change of electrical density of 
states (DOS) in low-dimensional structures,!^ such as semiconductor quantum wells, quantum 
wires, quantum dots, and single molecular devices.!^ Another approach is to search or design 
electron-crystal-phonon-glass materials by phonon engineering (phononics).!^ The idea is to re¬ 
duce the phonon thermal conductivity while keeping the electrical conductivity intact.!^ Recently, 
combining these two approaches, phonon engineering is used in low-dimensional structures, like 
silicon nanowires, to boost their ZT^^ A natural question then arises: Is it possible to do the 
opposite, utilizing low-dimensional electrical transport in bulk phonon-glass crystals? To show 
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this is indeed possible, we need to start from bulk materials with low thermal conduetivity. 

Thus far, the main focus of the thermoelectric research community is inorganic semiconduc¬ 
tor materials,including Bi 2 Te 3 ,^PbTe,^SiGe,^SnSe,^perovskites,^and Instead, 

we focus on organic semiconductors, especially molecular crystals, which have been intensively 
studied in the development of organic photovoltaic cells, ^light-emitting diodes,^ and field-effect 
transistorsCurrently, their thermoelectric properties have received more attention.Com¬ 
pared with inorganic materials, the advantages of organic thermoelectric materials are lightweight, 
flexible, cheap and easy to process over large areas Due to the weak bonding between different 
molecules, organic molecular crystals have natural low thermal conductivity. Recent experimen¬ 
tal progress has resulted in several orders of magnitude increase of ZT, from 0.001 qq^IMHI 
to 0.42.12^*22] Moreover, it is now possible to design and modulate the transport characteristics of 
organic structures.*^ Among them, Tt-Tt stacking molecular crystals have proved their superior 
electronic transport properties along the development of organic field-effect transistors.*^ Interest¬ 
ingly, it has been shown that, their electronic transport properties can be controlled by tuning the 
stacking angle and distance.*2^*281 

Here, we attend to Bis-Dithienothiophene molecular crystal (BDTMC) which has the interest¬ 
ing electrochemical and optical properties.*^ As shown in Fig.[^a), each dithienothiophene molec¬ 
ular possess three fused thiophene rings and they form crystal structure through k-k stacking. The 
way that small molecules stack together determines the structure and electronic property of the 
organic crystal. Herringbone stacking is one of the common stacking structures, whose electronic 
band structure is normally two-dimensional, as a result of n-n stacking extending in two dimen¬ 
sions. The BDTMC favors coplanar stacking*^ and has high mobility .*221421 Due to the strong Tt-Tt 
overlap only along one direction, it shows a quasi-one-dimensional (QID) band with a large band 
dispersion (Fig[2a-b)). This QID band structure is different from common two-dimensional band 
in molecular crystals and appropriate for our study on thermoelectrics. Besides, its thermoelectric 
transport properties are under experimental investigation.*^ 

Based on a general one-dimensional tight-binding model validated by Density Functional The- 
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ory (DFT), we employ the semi-elassieal Boltzmann transport theory to study its eleetronie trans¬ 
port properties. Meanwhile, the phonon thermal eonduetivity is ealeulated using the elassieal 
moleeular dynamies (MD) simulation. We show that, QID eleetronie transport and phonon-glass¬ 
like thermal transport are realized at the same time in bulk BDT moleeular erystal, whieh suggests 
a great thermoeleetrie potential for n-n staeking organie moleeular erystals. 

Results and Discussion 

Electronic structure. 

The lattiee and eleetronie strueture of BDTMC are ealeulated from Vienna Ab-initio Simulation 
Paekage (VASP).I^The Perdew-Burke-Emzerhof (PBE) version of generalized gradient approxi¬ 
mation (GGA) is used for the exehange-eorrelation potential.!^ The DET-D2 method of Grimme 
is used to take into aeeount the van der Waals (VDW) interaetions between different moleeules®! 
(details in Supporting information (SI), 1(A)). Eigure[2a) shows the relaxed strueture. It forms a 
trielinie Bravais lattiee, with two moleeules within one primitive unit eell, shown in dashed lines. 
The eonduetion (CB), valenee band (VB) strueture and DOS are shown in Eig. [^b) and (e). We 
observe a strong band dispersion along y direetion, G-B-G, eomparing to other direetions. This 
indieates a large overlap of n orbitals between moleeules along y direetion, due to the small angle 
between y and the n-n staeking (n) direetion. The band strueture along y is well re-produeed by a 
fitting using the dispersion of ID tight-binding model (red dots). Model details are given in SI, II. 

Due to a small eleetronie overlap between moleeules in the same unit eell, the eonduetion and 
valenee bands eome in pairs, eaeh mainly loealized in one moleeule. By eomparing the eharge 
distribution of the bands with the eleetronie states of the isolated moleeule pair, we find that the 
VB and CB are simply formed by the highest oeeupied moleeular orbitals (HOMOs) and lowest 
unoeeupied moleeular orbitals (EUMOs) of eaeh moleeule, respeetively [Eig.[2(d) and (e)]. The 
strong anisotropy in band dispersion suggests the formation of QID band strueture along y diree- 
tion. This is further supported by the van-Hove-like high DOS near the band edges [Eig.[2e)]. Eor 
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ID system, the DOS is inversely proportional to the hopping element t. Thus, the DOS at valenee 
band edge is around one-third of DOS at the eonduction band edge. Hereafter, we foeus on the 
eleetrieal properties along y direetion. 



Figure 1: The lattice and DFT band structure, (a) The relaxed strueture of the BDT erystal 
along different direetions. The lattiee veetors are: ai = (16.82, 0.087, 0.015) A, az = (16.40, 
3.85, 0.015) A, and as = (-2.08, -0.24, 10.68) A. The y-direction is approximately perpendieular to 
the moleeular plane, (b)-(e) DFT band strueture and density of states (DOS) of the BDT erystal. 
The energy zero is set to the top of the valenee band. The reeiproeal spaee eoordinates of high- 
symmetry points are G = (0, 0, 0), B = (-0.5, 0.5, 0), D = (0.5, 0.5, 0), and Z = (0, 0, 0.5), in unit of 
the reeiproeal lattiee eonstants. The red dots represent the bands fitted using one-dimensional band 
dispersion £ ~ —2tcos{kLo), with t = —0.2, —0.2 eV for the two valenee bands, and —0.05, —0.08 
eV for the two eonduetion bands (see text for details), (d) Charge density distribution of the two 
nearly degenerate HOMO orbitals in one isolated unit eell. They are loealized in eaeh moleeule, 
and form the two valenee bands in the lattiee strueture. (e) Corresponding LUMO orbitals forming 
the eonduetion bands. 


The intrinsie transport meehanism in organie moleeular erystals has been under aetive debate. 
Band and hopping models have been used to understand the experimental results.^^^l^It is eom- 
monly aeeepted that band-like model should be used if the inter-molecular hopping element t is 
larger or comparable to the molecule reorganization energy A .ESEH Following the method in Ref. 
I47l we obtained A of ~ 0.18 eV, comparable to t. Thus, we use the band model to study transport 
along y direction. This is also supported by the fact that hopping model predicts a mobility that is 
one-order-of-magnitude smaller than the experimental value in a similar molecular crystal.l^ 
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Based on the band structure obtained from VASP,I^ we use the BoltzTraP codeia to calculate 
its contribution to the transport coefficients, e.g., (Kg + S^oT)/x, S^a/x, with x the constant re¬ 
laxation time. Rigid band approximation is used to relate the Fermi level (Ef) position with the 
electron/hole doping concentration [Fig.j^a)]. The weak dependence of on doping within three 
decades is due to the van-Hove-like DOS at the band edges. Thus, similar to QID structure, we 
expect a high Seebeck coefficient shows up in the bulk molecular crystal. The DFT result is further 
compared with a ID model calculation using the parameters obtained from fitting the band struc¬ 
ture (Details in SI, 1(A)). We find excellent agreement between the two approaches near the band 
edges [Fig. I^b-c)]. This agreement validates the ID model, which we will use in the following. 

Electrical transport. 

A common approach to obtain the thermoelectric transport coefficients based on the DFT plus 
BoltzTraP calculation is to use a constant x. It is obtained either from fitting of the experiments, 
or from the theoretical estimation. To estimate the relaxation time, we need to consider different 
types of carrier scattering process. Defects, charge traps may scatter the carriers strongly, and 
even invalidate the band-like transport model used here. But these scatterings are extrinsic, and 
depend on the quality of the sample. The charged impurity scattering due to doping is another 
source of scattering. To take it account, we should consider the screening of the charged impurity, 
which is a problem that deserves separate study. Here, we only take into account the intrinsic 
scattering mechanism, namely the interaction of electrons with acoustic phonons, which is present 
independent on the quality of the sample. Due to the QID electronic structure, using the Bardeen- 
Shockley deformation potential theorywe go beyond the constant -t approach by taking into 
account the k dependence of x{k) (SI, II) 

1 ^. ksTD^ 
m ^ h^c\vk\' 
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Figure 2: DFT plus BoltzTraP calculation and model fitting, (a) The Fermi level {Ef) as a 
funetion of eleetron (red) and hole (blaek) doping at T = 300 K. The weak dependenee of the 
Ef on doping aeross three deeades is eharaeteristie of quasi-one-dimensional energy band, (b)-(e) 
The dependenee of K^jx = {Ki,-\-S^oT)/t and S^o/x on Ef. The points are results from DFT 
plus BoltzTraP ealeulation, while the lines are from the ID tight-binding model with the hopping 
element t fitted from the band structure in Fig. [^b). In (b), xyz means summing up all three 
directions, while y means contribution from y direction only. As can be seen, the contributions 
from X and z are negligible. 
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Here, kg is the Boltzmann eonstant, T is the temperature, C = d^E /(Lq^AL^) is the ID elastie 
eonstant, with AL = SL/Lq. D = Afi/AL is the deformation potential eonstant, with Afi the energy 
shift of the valenee or conduction band edge (SI, 1(B)). This makes our approach different from 
the common BoltzTraP calculation, and much closer to the real situation. The relevant transport 
coefficients Kg, o are shown in Fig. |^a) and (b). We see a better transport performance of the 
holes. It can be attributed to their larger hopping element t, leading to larger group velocity (SI, 
II). Moreover, a strong deviation of the Wiedemann-Franz law is observed in Fig. |^c), typical at 
the band edges or low-dimensional structures. The smaller value of Kg/{To) <L indicates reduced 
electron thermal conductivity, and good thermoelectric performance. 


Phonon transport. 

To get ZT, we also need the phonon thermal conductivity Kph, which was calculated by equilibrium 
MD simulation using LAMMPS^ (Details in SI, III). The following results are obtained from 
a simulation cell of 4 x 6 x 4 unit cells in x, y, and z directions, respectively, to overcome the 
finite size effect (SI, III). As shown in Fig. Q the thermal conductivity of BDTMC shows a very 
weak temperature dependence over the temperature range considered (100 - 350 K). This weak 
temperature dependence is normally observed at high temperature limit, well within the validity 
of classical MD simulation. We get a value of 0.34 ±0.02 W/m-K at 300 K, which falls into the 
common range of 0.1 ~ 1 W/m-K for organic molecular crystals 

The low thermal conductivity originates from the weak intermolecular bonding of VDW na¬ 
ture, in contrast to the strong intramolecular valence bonding. Different kinds of bonds result in 
very different frequency/time scales of inter- and intra-molecular dynamics. The thermal conduc¬ 
tivity of the molecular crystal is dominated by the low-frequency inter-molecular vibrations. Their 
frequency mismatch with the intra-molecular vibrations prevents efficient heat transport between 
molecules, good for thermoelectric performance. That is, most of the energy is localized inside 
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Figure 3: Thermoelectric transport coefficients at 300 K. (a)-(b) The electron thermal con¬ 
ductivity Kg and electrical conductivity <7 as a function of electron (n) and hole (p) doping, (c) 
Kg/{aT) as a function of doping, in unit of the Lorentz number L = 7r^kg/{3e^), such that the 
Wiedemann-Franz law corresponds to Ke/{oT) = L. (d) ZT as a function of doping. 
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each molecule, instead of transferred to other molecules. This can be seen from the large overshot 
of tCph as a function of correlation time in Fig. S4. 



Figure 4: Phonon thermal conductivity. The phonon thermal conductivity as a function of tem¬ 
perature obtained from the classical molecular dynamics simulation. 


Thermoelectric performance. 

We are now in the position of evaluating the thermoelectric figure of merit ZT. As is shown in 
Fig[^d), we get an optimal room-temperature ZT of 0.38 at electron doping of ~ 3.2 x 10^° cm~^, 
and 1.48 at hole doping of ~ 1.57 x 10^® cm^^. The corresponding Seebeck coefficients are —199 
/rV/K for electrons and 266 /iV/K for holes. The electron and hole mobility at optimal doping 
are ~ 2.4 and ~ 12.8 cm^/V-s, respectively. These results are summarized in Table The fol¬ 
lowing facts are note-worthy: (1) Although doping the organic semiconductors is still challenging 
in experiments, the optimal doping level we obtained here has been achieved in Pentacene experi¬ 
mentally;^ (2) A much higher carrier mobility has been experimentally observed in C8-BTBT;^ 
(3) Although the optimal electrical conductivity we obtained here has not been realized so far (the 
highest experimental value that can be achieved^ is ~ 11000 S/m), there is no fundamental dif¬ 
ficulty in realizing both (1) and (2) in the same material; (4) On the other hand, a recent study 
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shows that it is possible to reduce the optimal doping by increasing the carrier mobilityFinally, 
in Fig.|^ we plot the temperature dependence of ZT. It shows an approximately linear dependence 
on T. We have checked that ZT saturates at higher T (Fig. S6). Since band-like transport at high 
T is questionable, we did not show it here. The ZT value is promising over a wide temperature 
range (100 - 350 K). 

Table 1: Summary of important parameters and results for the conduction (CB) and valence band 
(VB): t the hopping element in the ID model, m* effective mass, D deformation potential, /r 
mobility at optimal doping, S Seebeck coefficient at optimal doping (n). 



t(eV) 

m* (mo) 

D(eV) 

n(10^'^cm ^) 

ZT 

At (cm^/V-s) S (AtV/K) 

Ke (W/m-K) 

CB 

0.05 (0.08) 

4.27 

-6.67 

3.19 

0.38 

2.4 

-199 

0.06 

VB 

0.2 

-1.5 

-8.55 

1.57 

1.48 

12.8 

266 

0.15 



Figure 5: Temperature dependence of ZT at optimal doping. 

The QID electrical transport and low phonon thermal conductivity we obtained here originate 
from the chemical nature of the inter-molecular bonding. In our calculation, the conduction and 
valence band structures of BDTMC are determined by only one parameter, t, reflecting the strength 
of the n-n VDW bonding. By comparing results of the CB and VB, we deduce that stronger 
bonding, thus larger t, results in better thermoelectric performance. 

Our results, together with recent development in chemical control of the stacking angle and dis- 
tance,!^^*^ point out possible way of further enhancing the thermoelectric performance by chemical 
or strain engineering of n-n overlap integral. DFT calculation confirms this idea. Within lattice 
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strain of 5%, we can achieve more than 20% change of t (Fig. S2). 

Before closing, we should mention that, as doping is increased, the impurity scattering rate is 
also increased proportionally to the dopant concentration, resulting in a decrease o and Kg. But 
their ratio should not change much. Considering ZT = 5 ^/{KejoT + Kpj^/oT), we expect the ZT 
to be smaller after including the impurity scattering. But to include it quantitatively, we need to 
consider the screening of charged impurity, which depends on the dimension and dispersion of the 
electronic band structure. Since the main idea of this work does not rely on these details, we leave 
them to future study. 

Conclusions 

To summarize, we propose a novel approach to search for highly-efficient thermoelectric materials. 
That is to explore or engineer low-dimensional electronic structure in phonon-glass bulk crystals. 
Through atomistic simulation, we have shown that the Tt-Tt stacking molecular crystals are par¬ 
ticularly suitable for such an approach. Both low phonon thermal conductivity and ID electronic 
structure originate from the n-n bonding. Besides proving the principle, our results also show the 
promising potential of 7l-7l stacking organic crystals as efficient thermoelectric material. Although 
we only considered organic crystals in this work, the idea should be equally applicable to inor¬ 
ganic materials. Indeed, during the preparation of this work we became aware of a recent work, 
exploring similar idea in inorganic compounds.^ 
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